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ABSTRACT 

We present SiFTO, a new empirical method for modeling Type la supernovae (SNe la) light curves 
by manipulating a spectral template. We make use of high-redshift SN data when training the model, 
allowing us to extend it bluer than rest frame U. This increases the utility of our high-redshift SN 
observations by allowing us to use more of the available data. We find that when the shape of the light 
curve is described using a stretch prescription, applying the same stretch at all wavelengths is not an 
adequate description. SiFTO therefore uses a generalization of stretch which applies different stretch 
factors as a function of both the wavelength of the observed filter and the stretch in the rest-frame B 
band. We compare SiFTO to other published light-curve models by applying them to the same set of 
SN photometry, and demonstrate that SiFTO and SALT2 perform better than the alternatives when 
judged by the scatter around the best fit luminosity distance relationship. We further demonstrate 
that when SiFTO and SALT2 are trained on the same data set the cosmological results agree. 
Subject headings: supernova:general 



1. INTRODUCTION 

The importance of modeling Type la super- 
nova (SN la) light curves is clear from the large 
number of methods that have been developed 
for this purpose. An incomplete sample in cludes 
MLCS/MLCS2k2 jgiess Press, fc Kirshnerl H99l 



Jha. Riess. fc Kirshnerl [20071) . stretch (iPerlmutter et al.l 
1997t IGoldhaber et alJ 120011 : IKnop et alj "" 
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stretch (IWang et al.ll2006|). Am is (lHamuv et aTJ[l 996b: 
Prieto. Rest, k Suntzefj 120061). B ATM (iTonrv et all 
20M) . CMAGIC (IWang et all (20031). SALT (|Guv et all 
20ull . and SALT2 (|Guv et all 120071 hereafter G07). 
For some SNe the different techniques produce quite 
different results, but when applied to a moderately sized 
sample the overall result s appear to be fairly similar 
(|Wood-Vasev et al.l 120071 ). These packa ges are most 
frequently used to derive relative distances which can 
then be used to constrain the cosmological parameters. 
Depending on the data set, current methods give relative 
distances to a precision of ~ 7 — 10%. The potential 
reward for any small improvement in how well distances 
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can be extracted from SN data is large. 

This paper presents the details of SiFTO, a light-curve 
analysis tool developed for use with data from th e Su- 
pernova Legacy Survey (SNLS; lAstier et al. 2006, here- 
after A06). SNLS follows a two-pronged analysis ap- 
proach: most of the steps are carried out twice using in- 
dependent methods, and the results are compared. This 
includes spectroscopic typing, photometry, calibration, 
light-curve fitting, and the extraction of the cosmologi- 
cal parameters. The details of the complementary light- 
curve fitti n g me thod (SALT/SALT2) are described in 
IGuv et alJ (pOOl ) and G07. By design, SiFTO takes a 
simpler algorithmic approach than SALT2, which should 
make the results somewhat more robust to the details of 
the training sample. SALT2, however, should have more 
growth potential as future data sets become available. 

The features of SiFTO are: 

• The shape of the light-curve is described by a single 
stretch parameter, the stretch in the rest-frame B 
band, sb- However, the model is generalized so 
that different observed filters stretch by different 
amounts as a function of wavelength and sb ( § 12 . 3j) . 



• SN la colors are handled by allowing the normaliza- 
tion of the fit in each observed filter to vary inde- 
pendently. The resulting colors are then combined 
to form a single color parameter when the results of 
the light-curve fit are transformed into a distance 

(§m. 

• For the A06 data sample, the RMS of the SiFTO- 
derived distances around the best fit cosmological 
parameters is 0.16 mag. This is identical to SALT2 
on the same data set, and is a clear improvement 
over SALT and MLCS2k2 (51731). 



• When we train SiFTO and SALT2 on the same 
data set using the same non-SN related settings 
(filter functions, etc.) the cosmological results are 
very similar (§[5]). 
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The outputs of SiFTO are a set of light-curve parame- 
ters (the time of maximum and the sb) and the spec- 
tral energy distribution (SED) of the SN as a function of 
time. From this information one can construct a distance 
estimate by finding the magnitude in the rest frame B 
band (tub) on a fiducial epoch and some set of rest-frame 
colors. The peak magnitude, stretch, and color can be 
combined to form a relative distance estimate. 

There are many similarities between SiFTO and 
SALT/SALT2. Unlike the other packages mentioned 
above, SiFTO, SALT, and SALT2 work in flux space 
by manipulating a model of the SED directly. By work- 
ing with the SEDs, the fits are naturally performed in 
the observer frame without the need for a distinct K- 
correction step. Like SALT2, SiFTO incorporates high- 
rcdshift SNLS data when deriving the model SED, which 
allows us to extend our model further towards the blue 
than is possible with only low-z data. This is particu- 
larly important when working with surveys that sample 
a large range of redshifts, since in this case some of the 
observed filters may be bluer than the rest-frame U band 
at the highest redshifts. This is true even of many sur- 
veys which obtain supplemental near-IR photometry. 

SNLS observes SNe in the range 0.2 < z < 1 in four 
filters (gMTMiMZ M), which are s imilar to the USNO sys- 
tem described in ISmith et alJ (|2002h . At the lowest- 
rcdshifts, these filters essentially sample BVR, and at 
highest %m i— > U and Zm l_ > B. ft is critical that our 
model handles both situations in a consistent fashion. 
For z > 0.55, gu is bluer in the rest frame than U, and 
the same holds for rjf at z > 0.9. Since the SN is of- 
ten still detected in these filters, if our model extends 
only to U (3600 A) we will have to ignore some of our 
signal. Incorporating vm for the highest redshift SNe is 
particularly helpful because the zm data are both nois- 
ier (due to fringing, sky brightness, and a drop in the 
CCD quantum efficiency) and harder to calibrate than 
9m, Vmi an d Im, and it is necessary that we have good 
quality observations of each SN in at least two filters to 
measure the color. 

There are many algorithmic differences between the fit- 
ters. SiFTO uses a spectral template as an input to the 
training process, although this template was informed 
by the results of initial SiFTO fits, while SALT2 in- 
cludes spectra more directly in the model derivation. The 
method of adjusting the SED models are also very differ- 
ent, with SALT2 in principle allowing for more detailed 
modeling of individual features. However, given the lim- 
itations of current training samples, the most important 
difference between SALT/SALT2 and SiFTO lies in how 
SN colors are handled. In SiFTO, the SED is adjusted 
to match the observed frame colors in each filter, rather 
than by imposing a single color model across all filters 
simultaneously as in SALT/SALT2. The justification is 
that SNe la show evidence for some additional variabil- 
ity in their colors, particularly in the near-UV, so forcing 
the fit to match a single color law at all wavelengths may 
give a poor result. Note that SiFTO, SALT, and SALT2 
do account for a relationship between stretch and color 
when the results are converted into a distance. In this 
paper, this relationship is taken to be linear. 

SiFTO is designed for use with modern, well measured 
SN data sets, and is not suited for the analysis of poorly 
sampled light curves. Specifically, if the number of obser- 



vations in a filter in the range -20 to 40 rest frame days 
relative to B maximum is less than about 3, then the 
fit in that filter may not be reliable. This proviso only 
applies to historical SN la samples - even at the highest 
redshifts, current observations are of high enough qual- 
ity that this is not a problem. Our approach is not well 
suited to the rest-frame near-IR (7 or redder), where a 
stretch-like prescription does not work well. 

Our base SEP model is an updated version of 
iHsiao et al.l ([20071 . hereafter H07), adjusted to match 
SN la observations in a number of synthetic filters (§ 13.2)1 . 
Our template includes an error snake representing the 
uncertainty in the template and intrinsic variability of 
SNe la (§ E3J). SiFTO is designed so that the effects 
of most systematics can be modeled by adjusting the 
color relation, and do not require fully retraining the en- 
tire template (§ 13.4)1 . The data used to train the SiFTO 
model is described in § [4j The derived light-curve param- 
eters are the most reliable when there are data within ~ 7 
rest frame days of the epoch of peak B luminosity (§E]). 

2. APPROACH 

2.1. Light- curve fitters and distance estimators 

The distinction between light-curve fitters and dis- 
tance estimators is often not made explicit in the litera- 
ture. The majority of the published packages are light- 
curve fitters, with the exception of MLCS/MLCS2k2 and 
BATM. A light-curve fitter attempts to find the best fit 
to a given set of SN la photometry. The parameters of 
this fit can then be converted into a distance estimate if 
desired, but technically this step is not part of the fit. A 
distance estimator attempts to find the distance directly 
rather than trying to obtain the best fit to the data. In 
both cases one is only trying to find relative distances; 
nothing described in this paper depends on knowing the 
absolute distance to any SN, the Hubble constant, or the 
absolute magnitude of any SN in any bandpass. 

The advantage of a distance estimator is that the out- 
put is directly what is desired for most applications of 
SN la data. Therefore, the products of such an analy- 
sis are simpler to use, and in principle such an approach 
may do a better job extracting the information directly 
relevant to this purpose. Their primary disadvantage is 
that, by their nature, they must use distance informa- 
tion in their training, usually in the form of residuals 
from the best-fit Hubble relation. This makes it diffi- 
cult to include both very low-redshift SNe (which are 
not in the smooth Hubble flow), and high-redshift data 
(where the residual depends on the cosmological param- 
eters). To use high-z data properly, it would be neces- 
sary to re-train the model from scratch for every value of 
the cosmological parameters one evaluates, which would 
be extremely computationally expensive. Therefore, in 
practice, light-curve fitters have access to a more data 
for training purposes. Neither approach is obviously su- 
perior, and we will not attempt to give a comprehensive 
comparison of their merits here. Because we find the ad- 
vantages of incorporating data from a range of redshifts 
in our training compelling, SiFTO is a light-curve fitter. 

2.2. SED modeling 

Like SALT/SALT2, SiFTO works by manipulating a 
SED model. Most fitters do not take this approach. In- 
stead, they are defined in some set of rest-frame filters 
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and use if-corrections to transfer the model to the ob- 
served frame. That is, they work in photometry space 
rather than SED space. This is simpler, but the fact 
that the if-corrections must be calculated and applied 
in a distinct step can lead to complications in ensuring 
that the two steps are fully consistent with each other. 
Working directly with the SEDs obviates this difficulty; 
formally, SiFTO, SALT, and SALT2 do not calculate K- 
corrections, although this is perhaps more of a pedagogi- 
cal distinction than a practical one. This choice requires 
that we work with fluxes rather than magnitudes, which 
is desirable anyways because using magnitudes will bias 
fits to low signal-to-noise data. 

We start with a base model for some fiducial SN la and 
then use additional parameters to adjust the SED to best 
match a particular SN. Two parameters are fairly obvi- 
ous (some overall flux normalization and overall time off- 
set), but beyond this the parameterization is somewhat 
arbitrary. Theoretical modeling of SNe la is sufficiently 
difficult that it offers little guidance here, so our choice 
of parameters is empirical. The existence of a relation 
between between the shape of the light curve and the 
peak brightness of the SN is very well established, as is 
the relation between the color and luminosity. These are 
in the sense that more slowly evolving SNe are observed 
to be more luminous (the wider-brighter relationship), 
and redder SNe are fainter (bluer-brighter). The for- 
mer is clearly related to SN physics, while the latter is 
very likely some combinatio n of SN properties and ex- 
ternal effects such as dust ([Branch fc Tammannl Il992t 
IConlev et ai1l2007f ). 

2.3. Light-curve shape 

The SiFTO model for light- curve shape is based 
on the stretch parame terization (|Perlmutter et ai]ll997t 
Goldhabc r et alJl2001f ). Here, the timescale over which 
a SN la evolves is modified by multiplying by a factor s. 
SNe with higher s evolve more slowly, with broader light 
curves. Changing the stretch simply widens or narrows 
the light curve in a linear fashion. 

In most previous applications of stretch either the same 
value was applied at all wavelengths or the stretch values 
were fit independently in each filter. The quality and 
quantity of SNLS data shows that the first is no longer 
sufficient, and we can improve on the second by relating 
the stretches of different filters. SiFTO uses an empirical 
relation between the stretch in B (sb) and the effective 
rest-frame wavelength of the observed filters. This is 
constrained to work like a simple stretch factor in rest- 
frame B. With this adjustment, the stretch model is 
a very good description of our observations, as will be 
demonstrated later. However, this approach does not 
work well in the rest-frame / and redder, which may be 
an issue for some data sets, although not SNLS. The 
SALT model is based on the stretch prescription in B, 
but not in other filters, while SALT2 allows for more 
general behavior. 

2.4. SN colors 

Rather than imposing a relation between different ob- 
served filters during the light-curve fit, in SiFTO this 
relation is used when the fit parameters are converted 
into a distance. The normalization of the fit to each ob- 
served filter is allowed to float independently during the 



fit. The SED is then adjusted using splines to match 
these observed colors. We find that single parameter 
functions (such as dust laws or the SALT/SALT2 color 
relation) are generally unable to reproduce the full va- 
riety of observed colors across more than two observed 
filters, which is likely indicative of some intrinsic vari- 
ability in SNe la. 

Therefore, the parameters in each fit are the time of 
maximum luminosity in the rest-frame £?-band to, the 
stretch sb, and Nf normalization parameters, one for 
each observed filter whose rest frame wavelength lies 
within the range our model is defined. Nf = 4 for most 
SNLS data, except at the redshift extremes. The results 
can be used to produce a corrected distance which makes 
a better standard candle: 

m corr = m B + a (s B - 1) - f3C. (1) 

a and (3 must be determined empirically along with the 
cosmological parameters, raj is the peak rest-frame mag- 
nitude in the B band, and C is some sort of color param- 
eter. In principle we could have several color parame- 
ters, one for each observed color, but it would be difficult 
to ensure that the color was being applied consistently 
across the entire sample. We therefore combine the dif- 
ferent observed colors into a single parameter, using the 
SN data themselves to derive the relationship between 
different rest-frame colors. 

One of the consequences of this relation is that any lin- 
ear variation of the intrinsic color of SNe la with stretch 
is absorbed by the a term. Therefore, our method does 
model some aspects of intrinsic color and how it relates 
to light-curve shape. The differences between this ap- 
proach and that currently used by MLCS2k2 are three- 
fold: First, MLCS2k2 includes a quadratic term in the 
light-curve shape. This has a small effect except at the 
extremes of the SN population. Second, MLCS2k2 makes 
use of several priors on the parameters, of which the color 
prior is the most important. Third, it assumes that any 
residual color that remains after the correction for light- 
curve shape is due to Milky Way-like dust, and there- 
fore j3 is (usually) set to 4.1 for the B — V color. In 
effect, these assumptions mean that SN colors can only 
derive from three sources: the shape of the light curve, 
dust, and possibly some scatter in the intrinsic colors 
which has no effect on the luminosity. MLCS2k2 could be 
modified to relax this set of assumptions. Therefore, the 
current distinction between intrinsic and extrinsic colors 
in MCLS2k2 versus SiFTO and SALT/SALT2 is more 
practical than philosophical, and is mostly related to the 
exact nature of the color and light-curve shape priors 
used by MLCS2k2. Whether or not these priors are ben- 
eficial depends on whether the smaller statistical errors 
obtained by using them outweigh the enhanced sensi- 
tivity to evolution or other systematic errors in the prior 
distributions (e.g., the distribution of extinction and how 
it interacts with survey selection effects) that they intro- 
duce; 

IConlev et al.l (|2007f l demonstrates that either dust 
along the line of sight to local SNe is either quite dif- 
ferent than we expect, or that SN intrinsic colors are 
not predicted completely by a linear or quadratic light- 
curve shape model by showing that (3 is significantly less 
than 4.1. We have further investigated this question by 
introducing a quadratic stretch term into our corrected 
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magnitudes, and find that it has virtually no effect on [3. 
In other words, adding an additional quadratic term like 
the one used in MLCS2k2 does little to address the issue 
of intrinsic versus extrinsic color. 

We build a model that can combine the different mea- 
sured colors for individual SNe into a single color parame- 
ter by taking SNe with simultaneous rest-frame UBV ob- 
servations and using this data to predict B — V as a func- 
tion of U — B and sb- Here, as elsewhere in this paper, 
colors are always measured at the epoch of B maximum. 
In effect, this is loosely similar to the SALT/SALT2 
approach, except that the color model has been sepa- 
rated from the light-curve fit and the SED model build- 
ing. In A06 a similar relation, implicitly included in the 
SALT training process, was determined purely by low- 
z SN data. U dat a is notoriously difficult to calibrate 
(e.g.. lBessell| [l990). and relatively little is available, al- 
though the situation has impr oved considerably recently 
(jJha. Riess. fc Kirshner|[2007f) . By using SNLS data we 
can overcome this problem - the rest-frame U generally 
maps to either observer frame gu or tm, both of which 
are well calibrated. In addition, this data is available for 
virtually all of our SNe. As a result, the color relation 
is considerably better constrained than was the case in 
A06. 

Our color relation takes the form 

(B - V) v _ Bs =a(U-B + 0.5) + b (s B - 1) + c. (2) 

We then combine this with the actual measured B — V 
(if available) using a weighted average to form our color 
parameter C, which therefore represents the B — V color 
of the SN at peak, usually after correction for Milky- 
Way extinction. This differs by an arbitrary constant 
offset from the definition used in SALT/SALT2. 

The data suggests that there is some additional intrin- 
sic dispersion in how well B — V can be predicted by 
U — B beyond that attributable to measurement errors 

in other words, the different colors do carry some in- 
dependent information. This should not be entirely sur- 
prising, as there has long been a theoretical expectation 
that SNe la are not as homogeneous in the near-UV as at 
longer wavelengt hs (iLentz et al.ll2 000) . This is supported 
by observations (|Ellis et al.l I2007T ) ■ We account for this 
by degrading the assigned uncertainty in the predicted 
B — V by an additional factor er c i r added in quadrature 
with the measurement error. This plays a similar role 
to the if-correction uncertainties discussed in G07, al- 
though it is expressed in a different fashion. 

Because of this additional scatter it is important that 
the statistic used to measure a, b, c allows for dispersion 
beyond that represented by the meas urement err ors. We 
use the multiple regression method of lKellvi (|2007f l. which 
also provides estimates of <7 c i r , and couple these fits with 
iterative 2.5a outlier rejection when deriving the relation. 

A further refinement is to extend this model bluer than 
the rest- frame /7-band. We define an art ificial filter de- 
fined as the £/-band filter of lBessell (jl990[ ) blueshifted by 
z = —0.2, which we refer to as Uq2, and which has an 
effective wavelength of ~ 3000 A. We then use a similar 
relation to the above to predict B — V given U02 — B. Not 
surprisingly, this relationship is not as well determined 
since it is measured with higher-redshift, and hence nois- 
ier, observations. Nonetheless, the additional informa- 
tion provided by t/02 is useful with the highest-redshift 



SNLS observations (fj5]). We denote the coefficients of 
this relation as 002 &02 and C02 , with a corresponding 
co2ch - The same technique can also be applied to V — R. 
Since this is not useful for the majority of SNLS data 
(which do not have rest frame R observations), we do 
not discuss these coefficients here. 

This gives us (potentially) three estimates of B — V - 
that from U — B, Uq2 — B, and finally the actual B — V. 
We form a weighted mean of whatever subset are avail- 
able to form C, taking care to propagate the covariances. 
The additional scatter represented by cr c i r and similar 
terms usually has a significant effect on the final weights. 
Large extrapolations are undesirable in this process be- 
cause they increase the importance of the boundary con- 
ditions used when we adjust the SED. To avoid this, we 
require that each of our synthetic rest frame magnitudes 
(U02U BV) has an observed frame filter within 650 A (in 
the rest frame), and do not allow the same observed fil- 
ter to be matched to multiple rest-frame filters. With 
the requirement in place, our fits are very insensitive to 
the form of the boundary conditions. 

The model blue-ward of U (~ 3600 A) is completely 
determined by SNLS data. This means that the question 
of how well this part of the model is calibrated observa- 
tionally is determined by our and tm calibrations, 
which are secure. The translation from magnitudes to 
flux is not as precise in this wavelength regime (i.e., we 
do not know the flux of Vega in C/02 as well as in the stan- 
dard filters) , and so in an absolute sense our model SED 
has some additional calibration uncertainties here. How- 
ever, because we do not have any observer frame data in 
this wavelength range, these are completely irrelevant to 
our current application. That is, the SED of Vega could 
be replaced with any function below 3000 A and it would 
have no effect on our derived distances. Of more interest 
is the concern that SNe la may not be good standard 
candles in the near-UV, and that there may be external 
factors (such as dust) which vary more with environment 
at 3000 A than at 3600 A. 

3. DETAILS OF THE MODEL 

Given our SED model <j> (described in § 13 . 2[) . param- 
eterized by the observed epoch t, the stretch in the B 
band sb, the epoch of peak luminosity in the B band 
to, the Nf phase- independent multiplicative factors rij, 
we define the effective broadband integrated flux Fj in a 
filter j of a SN observed at a heliocentric redshift Zhei by 

JiT ( idir7> *' *0' n i> SB ) t j ( a ) AdA 

p. = v " c ' L (i\ 

(l + z hel )J °°Tj(\)\dX ' [> 

where Tj(X) is the dimcnsionless total (atmosphere + op- 
tics + filter + CCD) system transmission function. Note 
that this is not the energy flux absorbed by the detec- 
tor, differing both in normalization and by including an 
extra factor of A inside the integral. The utility of this 
quantity stems from its relation to astronomical magni- 
tudes: m = — 2.51og 10 Fj + const. The factor of A is a 
result of the fact that astronomical magnitudes are de- 
fined in terms of counts, not energy. The calibration 
process by which observations are placed on a standard 
photometric system involves a similar term over some 
reference spectrum whose magnitudes are known; this is 
conventionally taken to be the spectrum of Vega. We 
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require the rij to be positive. The Tj are never known 
perfectly, and this should be included in the systematics 
error budget of any cosmological analysis, including the 
consequences for SiFTO or any other light-curve fitter. 
An exploration of how well these are known for differ- 
ent supernova samples is beyond the scope of this paper. 
We expect to update SiFTO as our understanding of var- 
ious calibrations improves and new SN samples become 
available. 
The light-curve fit then minimizes 



2 _ K- t j\ t > to, ni, SB, Zhel) - Jij) 



j=i i=i afj + <y 2 F it, to, rii, sb, Zhei, X 3 cS ) 

th 



(4) 



where Nj is the number of datapoints in the j filter, o~p 

is the template error, X J eS is the effective wavelength of 
filter j (see below), and fa are the observed data points 
with errors o~ij . 

The parameterization of (f> is the most critical element. 
We start with a model for the SED of a fiducial sb = 1 
SN la, 0o (A, t), where we define sb so that the mean 
value averaged over our training sample is unity, and 
0o is a time series of spectral templates covering the 
range 2000-9000 A with a unique SED at every epoch 
from —19 to +70 days. We then form the effective epoch 
r = (t — to) / (sb (1 + ^hei)), an d then use this quantity 
to select the appropriate SED from <po, linearly interpo- 
lating between the two nearest tabulated epochs. 

We next apply the stretch model as described below, 
and adjust the SED to match the observed colors so that 
the final result is a single SED which can be used to pre- 
dict the magnitudes in any filter at any desired redshift. 
Usually this is used to compute rest-frame magnitudes 
in some standard set of filters such as U02UBV. 

The method used to adjust the SED is similar to that 
described in H07: we multiply the SED by an interpolat- 
ing cubic spline under tension (<S), with the spline knots 
placed at the effective wavelengths of the observed fil- 
ters. In order to ensure that the SED remains positive, 
we work in log space. An example is shown in figure [TJ 
We experimented with a variety of functional forms to 
smoothly scale the spectrum in place of splines, and 
found similar results using linear or higher-order poly- 
nomial interpolation. 

3.1. The stretch model 

In the SiFTO model, Sj, the stretch in filter j, is not 
necessarily the same as sb ■ Our model for this correction 

cff 1 , where A^ and X J eS are the 

effective rest frame wavelengths of B and the j th filter 
for our fiducial, Sb = 1 SED at t = t$: 



is Sj = h (s B , Af ff - X{ 



A cff 



J0o {\t = to)T B {\) A 2 dA 
J Tb (A) AdA 



and similarly for X J eS , although for the latter the redshift 
of the SN needs to be taken into account. Given the 
stretch in the observed filter, we adjust the corresponding 
flux normalization by multiplying it by the ratio: 



TZj (t, to, SB, Sj, Zhel) 



AdA 



S B (l + 2hel) 



Tj (A) AdA 



This gives observer-frame flux in each observed filter. 
There is still only one light-curve shape parameter, sb- 

The stretch function h is represented by a spline, and is 
derived from both external and SNLS SN la data as fol- 
lows. We require that at least one of the filters lie within 
250A of rest frame B in order to ensure that sb is well 
determined. However, in addition to fits to individual 
filters, we also consider nearest neighbor combinations of 
the SNLS filters (qm +Tm, t m and %m + zm), and 

include the SN if the mean wavelength of these filters 
meets the above requirement. In these cases B is well 
bracketed, so sb is also well constrained. The stretch 
and (mean) rest-frame effective wavelength is recorded 
for each fit. The results are shown in figured] 

Alhough the relationship between sb and Sj varies 
slightly with wavelength, a general trend is apparent: 
sb is typically smaller than Sj for small sb and larger 
for large sb ■ It is also apparent that the low-z and SNLS 
samples probe different ranges in sb (with the low-z sam- 
ple containing more low-s^ SNe and the SNLS sample 
more high-s^ SNe) . This is not surprising given expected 
selection eff ects and the possib ility of an evolving popu- 
lation (e . g ., EoWriEaIJ[l003) • 

To model the trends, we fit the relation between sb 
and Sj using least-squares in ~ 100 — 200A wide bins, us- 
ing a cubic spline with three knot points in sb smoothly 
joined to a linear relation at low and high sb- Each 
wavelength bin contains 20-60 datapoints; three exam- 
ple bins corresponding to U, V and R are shown in fig- 
ure [H The relationship between sb and Sj at any inter- 
mediate Ag ff can be found by linearly interpolating the 
fit coefficients in wavelength. We do not probe all areas 
of A e ff equally well, with clusters around the low-z filter 
mean- wave lengths; this is an area where larger samples 
of SNe la will be extremely beneficial. Examples of the 
derived function are shown in figure 0J 

We can test the necessity of this more complex treat- 
ment of stretch by applying the SiFTO model to the 
photometry from A06 (see §7.31 for more details) . For 
the model described above, we find a \ 2 of 2890 for 2822 
degrees of freedom. If we remove the wavelength depen- 
dence of the stretch model the \ 2 climbs to 3037, and 
if we instead use a purely linear model in place of the 
splines, the \ 2 is 2935. These correspond to a decrease 
in likelihood by a factor of 70 and 3, respectively. How- 
ever, we note that the cosmological results and the scat- 
ter around the best fit are almost identical in all three 
cases ( < 0.001 mag difference in the RMS). The lack of 
improvement in the relative distances can be interpreted 
as an illustration of the difference between a light-curve 
fitter and a distance estimator. 

3.2. Constructing the SED model 

The starting point for our base SED model 0o is an 
updated version of the time series of spectral templates 
of H07 which incorp orates the addit i onal low-z spec- 
troscopy presented in iMatheson et alj ([2008) . The pro- 
cedures of H07 are designed so that the relative strengths 
of adjacent features are correct, but are not intended to 
adjust the broadband colors or their evolution with time, 
so some modification is necessary. We developed a set of 
light-curve templates that are used to address this issue. 
These templates are purely an intermediate step, and 
play no role in the fitting process once 0o is determined. 
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There are many approaches one might use to build 
the templates. For SiFTO, we construct a set of syn- 
thetic filters and then use SN data to build light-curve 
templates in each of these filters. We then adjust <fio to 
match these templates using interpolating splines as de- 
scribed earlier. This process requires ^-correcting the 
observed data points, which are at a variety of redshifts, 
to the nearest template filter. Since the A"-correction de- 
pends on the model SED, this procedure is carried out 
iteratively. We follow a fairly standard approach by pa- 
rameterizing each light-curve template as a set of cubic 
splines with the kn ot points tuned by hand to best de- 
scribe the data (e.g., [Goidh^eriESJ(l2Q0i|)). The num- 
ber of knots ranges from 3 (for poorly constrained filters) 
to 8 for the redder filters which have a more complex tem- 
plate. Since cubic splines are used, the templates are only 
smooth up to first order in the derivatives; applications 
which require higher degrees of smoothness should not 
use SiFTO. The splines are joined smoothly to a polyno- 
mial at early times, and an exponential decay function at 
late times. The location of the join can be set for each fil- 
ter; generally the early time polynomial joins the splines 
around epoch —10, and the exponential tail around +50, 
although in some filters this is much lower due to the lack 
of good late-time data. The early-time polynomial is of 
the form / = a p (t- t cxp ) 2 +b p (t- t cxp ) 3 , where t cxp is 
the date of zero flux (i.e., the date of explosion). There- 
fore, for our template, the colors of the SN at early times 
are not constant in all cases due to the cubic term. We re- 
quire the polynomial coefficients to be positive to ensure 
that the early time light-curve increases monotonically at 
early times, which we expect on physical grounds. Linear 
and quartic terms were not found to improve the fits. 

The choice of which filters to construct templates in 
is fairly arbitrary. If only low-z data is used for train- 
ing, then a natural choice would be the Landolt U BVRI 
filters, redshifted slightly to the mean redshift of the sam- 
ple. Since we use both low-z and SNLS data, this is no 
longer optimal. We want to make use of information 
blue-ward of U, and it is also useful to insert filters be- 
tween the more widely spaced Landolt filters to make 
better use of our observations and improve the smooth- 
ness of the model. 

We use a set of synthetic logarithmically spaced filters 
between 4800 A and 6800 A, and also include B and 
/. In addition, we include three bluer filters which are 
blucshifted and widened versions of the first log filter. 
The B filter is included to make it easier to compare the 
resulting template with previous work, and I to ensure 
well-behaved boundary conditions while building 6q. We 
do not use the part of the model redder than 7100 A when 
applying SiFTO to observations because our stretch pa- 
rameterization is not expected to work well. The syn- 
thetic log filters are denoted as LOG5 - LOG8, and the 
three bluer filters are ARB0 - ARB2. The light-curves in 
ARB0 and ARB1 are determined purely by SNLS data. 
Our filters are shown in relation to the Landolt filters in 
figure[5l ARB2 matches U (but is slightly redder), LOG6 
roughly matches V, R is split up into two narrower fil- 
ters which are better matches to the SNLS observations 
(LOG7 and LOG8), and there is an additional filter be- 
tween B and V (LOG5). The filters are summarized in 
table [TJ The SiFTO model should only be used for ob- 



servations from 2700 to 7100 A (rest frame). 

We use cubic terms in the early-time polynomials for 
the ARB0, ARB2, and LOG5 filters. For ARB0, includ- 
ing this term reduces the \ 2 °f the template fit to our 
training sample from 171 for 155 degrees of freedom to 
167. For ARB 2 the \ 2 goes from 4613 to 4515 for 1605 
degrees of freedom, and for LOG5 from 2280 for 1004 
degrees of freedom to 2001. We caution that these % 2 
values do not incorporate the template error snake as 
described in H3.31 and hence we do not expect the re- 
duced x 2 to be close to unity. It is difficult to evaluate 
how significant the improvements in the \ 2 are, since the 
process of building the error snake involves inflating the 
errors to take into account the residuals from the model. 
However, in the other filters the best fit value of the cubic 
coefficient is zero. 

The procedure of building the templates in each filter 
is as follows: First, all of the data is fit using SiFTO with 
the current version of (j>o- The observed data is K- and 
MW-extinction corrected to the best matching template 
filter, and the timescale is normalized using to and the 
stretch in each filter Sj. The data points from each SN 
in each filter are normalized to an arbitrary value in or- 
der to remove the colors of the individual SNe from the 
template. The templates are then fit to the datapoints in 
each filter, and O is adjusted using the spline procedure 
to match the model. The procedure is then iterated un- 
til convergence, which typically requires 3 — 5 iterations. 
We also adjust the template so that t = to represents the 
peak flux in B, and adjust sb so that the mean over our 
training sample is unity. 

It is convenient for our template to have roughly the 
colors of a fiducial SN la. Given the asymmetry of the 
color distribution we use the median color to define these 
values. There is no single SN with exactly this color, and 
if there were it would not constrain all of our template 
filters. Therefore, we take the mean 10% of the SN closest 
to the median color to define each color, and use this in 
the normalization process above. Note that this process 
does not affect the resulting light-curve fits as long as 
sharp, unphysical gradients are avoided. Changing the 
colors individually by up to 0.3 mag has no effect on 
the resulting fits because our procedure adjusts the SED 
template to match the observed colors. 

With the above ingredients, we can write cf> as 

(j)(X,t,t ,ni,s B ) = </>o (A, — t j— — — r ) 
V s B (l + z ho i)/ 

x S (X;TZi (i,io,SB,Sj,-Zhci)^) (5) 

where S is the spline used to make the SED have the 
desired colors, and not that used to construct O . S 
is uniquely determined by the combination of the TZiUi 
products, cf>o, and the observed filters. <j) can then be 
used to evaluate the flux in any desired filter, and not 
just those actually observed. 

3.3. The template error snake 

SiFTO incorporates an estimate of the uncertainty in 
<Pq. This represents both the uncertainty in our modeling 
process, and, much more importantly, some measure of 
the intrinsic variability of SNe la. SALT2 and MLCS2k2 
also incorporate an error snake, but some older fitters do 
not. The uncertainty in the template is very important 
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TABLE 1 
Template Spline Parameters 



Filter 


w 


A c ff a 


"T~>T t 7TTH It 

1 WHM 


. b 
1-poly 


N s c 


Hail 


ARBO 


2689 


2798 


529 


-11.5 


3 


12 


ARB1 


3187 


3206 


530 


-11 


4 


25 


ARB2 


3685 


3704 


510 


-11 


5 


30 


B 


4413 


4330 


893 


-10 


6 


50 


LOG5 


4784 


4755 


497 


-10.5 


5 


41 


LOG6 


5382 


5363 


559 


-10 


6 


48 


LOG7 


6055 


6030 


628 


-10 


6 


50 


LOG8 


6812 


6854 


707 


-11 


7 


47 


I 


8060 


8008 


1526 


-11 


7 


17 



Note. — Information about the template niters and the 
functional form used to fit each filter when deriving <po- Ex- 
cept for B and 7, these niters are synthetic. These values are 
essentially hand-tuned. 

a The effective wavelength for this filter when observing a fidu- 
cial SN la at z = 0. b The epoch of the join between the splines 
and the polynomial early form, relative to the epoch of B 
maximum. The number of internal knots in the spline between 
*poly an d Hail d The epoch of the join between the splines and 
the exponential late form, relative to the epoch of B maximum. 

when fitting low-z data, where the photometric uncer- 
tainties can often be quite small. It also plays a critical 
role when evaluating gM and tm observations of high-z 
SN, where the model is more poorly constrained. We 
tabulate the template error at the wavelengths of each 
of the template filters in one day bins and then linearly 
interpolate between these errors in wavelength space to 
get the error in the SED, which gives us the of term in 
equation 2] 

There are several methods used to derive the error 
snake, which give similar results. The most direct ap- 
proach is to fit all of the SNe using SiFTO, and then 
determine the additional error that would make the % 2 
of the data relative to the model be 1 per degree of free- 
dom. The difficulty with this approach is that at some 
epochs the measurement errors of individual observations 
dominate the scatter. The SNLS error bars are roughly 
constant with epoch in flux space because it is a rolling 
search, but a 10% variation in the template flux at any 
given epoch is much larger in absolute flux at peak than 
at late or early times. Therefore, measuring the same rel- 
ative variation at late times is more difficult, requiring 
an accurate cancellation of two numbers subtracted in 
quadrature. Compensating for this fact is the tendency 
of the template error to be larger in relative terms at late 
times which partially arises from the fact that we have 
renormalized the fits to each SN at peak. In order to im- 
prove our sensitivity at late and early epochs, we adopt a 
slightly modified procedure. In the core of the light curve 
(—10 to +15 days) we use the x 2 method directly, but 
at late and early times we substitute a measure based 
on the RMS, but multiplied by the ratio between the x 2 
result and the RMS as measured in the core. The results 
of this procedure are generally quite similar to that of 
the raw \ 2 method (within about 10%). 

The other method we consider is to use bootstrap- with- 
replacement to estimate the template errors. This gives 
similar results, but the x 2 + RMS approach is more ro- 
bust and slightly more conservative in that it gives larger 
errors (by ~ 20% on average), so we use this for our final 
model. The error snake is used iteratively using the tem- 



plate construction process. Some examples are shown in 
figure [6] 

Since we usually have a large number of observations 
in each epoch bin for each filter, the formal statistical 
error in how well we can measure the mean template is 
much smaller (~ 100 times) than the template error we 
quote. Some of this is probably due to failures in the 
photometry, uncertainties in the if-corrections, and in- 
consistencies between different samples. However, look- 
ing at any of these in detail (for example, by considering 
only data from a single source such as SNLS, or looking 
at SN in a narrow redshift range where if-correction un- 
certainties are less important in a relative sense), none 
of these effects is the right order of magnitude to explain 
the observed scatter. We conclude that SN variability is 
the dominant factor in the template error, except in our 
bluest artificial filter where the sample is small. 

Like SALT2, our error snake does not incorporate off 
diagonal terms - i.e., the residual from the model at a 
given epoch is assumed to be uncorrelated with nearby 
epochs. Given the complexity of our model space, the 
current data set does not provide robust constraints on 
these terms. MLCS2k2, which has a simpler model and 
does not work with a SED, includes off-diagonal terms 
with simplifying assumptions. 

3.4. Modeling systematics using SiFTO 

Modeling of systematic errors is becoming increasingly 
important for SN cosmology. One issue which has been 
neglected is that because our SN model is empirically de- 
rived from SN data, the effects of the systematics on the 
model must also be considered in the error budget. This 
is important even if no high-z data is used in the train- 
ing. An advantage that SiFTO has over most other light- 
curve fitters is that the effects of epoch-independent sys- 
tematics can be studied simply by re-deriving the color 
relation, rather than having to fully retrain all aspects 
of the model. For example, a change in the SNLS gM 
zeropoint data does not affect the SED model because 
during the training process the data from each SN is 
normalized to an arbitrary level. It will, however, affect 
the relation between the measured colors of the SN popu- 
lation. Approaches which fold the modeling of SN colors 
into the model derivation more directly (SALT/SALT2, 
MLCS2k2) do not have this luxury. 

This does not apply to any effect which can change the 
shape of the light curve, such as uncertainties in the filter 
responses. However, in most cases these are second order 
effects compared with the changes in the derived color 
relation. This makes systematic analysis with SiFTO 
comparatively easy. 

4. DATA 

Our training sample consists of photometry for nearby 
SNe la from the literature as well as high-redshift SNLS 
observations. The low-z data generally has higher signal- 
to-noise ratios, but are also much more heterogeneous. 
More importantly, the SNLS data offers considerably 
better wavelength coverage than the low-z sample, both 
in the near-UV and between the standard Landolt filters. 
Furthermore, the rest-frame U band data from SNLS is 
better calibrated and more reliable than that from the 
low-z data. 
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The photometric references of the 72 low-z SNe used 
in this paper are given in table O We have not included 
sub-luminous (SN 1991bg-like) SNe la in our sample be- 
cause they are not well represented by our base SED, and 
have excluded extremely peculiar SNe such as SN 2002cx 
(|Li et all 12003). Furthermore, we have not included all 
of the available photometry for every SN. In cases where 
photometry is available from multiple sources, we com- 
pare the light-curves from multiple bands and, if they are 
not consistent, choose some subset of the sources. We try 
to choose the photometry that has the best coverage, or 
that from the data source with the most other SNe (i.e., 
we pre fer observations from a large sample like lJha et al.l 
(2006)). We also exclude individual bands on a SN by 
SN basis where the photometry seems to be internally in- 
consistent or has additional problems (e.g., the U band 
data for SN 1999ee). 

In addition we use 98 high-redshift SNLS SNe from 
the first three years of the survey. Our training process 
requires that the peak of each observed filter be fairly 
well constrained so that we can re-normalize the light- 
curves accurately, a necessity since we do not impose a 
color relation during the training process or fits. There- 
fore, to use an individual bandpass for a SN, we require 
that it have at least one datapoint within 5 rest-frame 
days of maximum. This requirement is far more stringent 
than is necessary to carry out a cosmological analysis, so 
less than half of our high-z data are used in the training 
process. We use only SNe with firm spectroscopic iden- 
tification. The H07 template makes use of SNLS spec- 
tra, which are particularly critical when constraining the 
near-UV SED. As was the case for the low-z sample, we 
exclude SNe which are known to be unusual, such as th e 
super-Chandrasekhar mass 03D3bb (iHowell et all feOOfA 

For SNLS, we adopt the radially-dependent SAGEM 
filter scans for the MegaCam filters, corrected for the f/4 
converging beam and incorporating the Mauna Kea ex- 
tinction curve at the mean survey airmass of 1.2. A06 
used radially averaged versions of the se scans. For the 
low-z data we adopt the iBesselll (| 1990T ) realization of the 
Landolt filter functions. Ideally we would have natural 
system magnitudes and bandpasses for all of the low-z 
photometry, but this information is generally not avail- 
able. Hopefully this will not be the case for future low-z 
data sets. Currently we do not implement any airmass 
dependence in our assumed filter functions, since tests 
show that this is negligible in comparison with the other 
bandpass uncertainties. 

5. RESULTS OF THE TRAINING 

The light-curve templates can be seen in figure [7] The 
different manner in which the SNLS and low-z data sets 
sample wavelength space is clear - only the SNLS data 
offer any constraints blue- ward of U (ARBO and ARB1), 
and the wavelength region between B and V (LOG5) and 
V and R (LOG7), while the nearby data provide most 
of the constraints in RI. However, in the filters where 
substantial coverage is provided by both sets (LOG3, B, 
LOG6), the agreement is impressive. In addition, the 
homogeneity of SNe la light curves is also apparent. Note 
that in the / band extensive culling has been performed 
because many SN do not match the template. This is 
acceptable because we are only including I to ensure well 
behaved boundary conditions, and do not use this part 



TABLE 2 
Low-Redshift SNe Ia Used to 
Build the Rest-Frame SED 
Model 



Name 


Filters 


Source a 


1981D 


UBV 


H91 


1986G 


UBV 


P87 


1989B 


UBV 


W94 


1 QQON 


UBV 


L98 


1990O 


BV 


H96 


1 QQO^f 


BV 


H96 


1 992ae 


BV 


H96 


1 992al 


BV RI 


H96 




BV RI 


H96 




BVI 


H96 


1 QQ9hh 


BV I 


1196 


1992M 


BV 


H96 


1 992ho 

J- ij £1 Vj \J 


BV RI 


H96 


1 QQ9hr> 


BVI 


H96 


1992br 


BV 


H96 




BV 


H96 


1993B 


BV 


H96 


1993H 


BV RI 


HQfi A 04 


19930 


BV I 


H96 




BV 


H96 


1 QQ4A/T 


BV RI 


R QQ A 04 


1994S 


BVI 


R99 




BV RI 


R05 A 04 


1995D 


BV RI 


R99 


1995E 


BV RI 


R99 


1995ac 


BVRI 


R99 




BV RI 


R99 




BV RI 


R 99 A 04 

-LLC/;?} nui 


1996C 


BV RI 


R99 


1996X 


UBVRI 


RQQ SOI 


1 QQfiah 

J_ U U U el 1 J 


BV 


R99 


1 99fihl 


BV RI 


R99 




BV RI 


R99 


1997E 


UBVRI 


J 06 




UBVRI 


J 06 


1997dg 


UBVRI 


J 06 


A-iJiJ l U-U 


UBVR 


J 06 


1998V 


UBVRI 


J 06 




UBVRI 


R05 




UBVRI 


TQQ S99 




UBV 


J 06 




UBVRI 


J 06 




UBVRI 


106 A 04 K00 

JUUj Au^t. IVUU 


1 QQQar 


UBVRI 


J 06 


1 QQQaw 

-L £7 C/ iyaVV 


BV RI 


S02a 




UBVR 


lOfi Kflfi 


1 QQQrl 


BV 


lOfi KTlfi 


1 999Hk 


UBV I 


KOI A 04 


1 QQQrln 


UBVRI 


J 06 


1 QQQpp 


BV RI 


S02b 


1 QQQpU 


BV RI 


lOfi K04 


lyyygu 


UBVR 


J 06 


1 QQQpti 


UBV I 


106 KOI 




UBVRI 


K04 


9nnnm 

^KJxJKJ^ It 


UBVRI 


J 06 


2nnnHk 

i-UwlK 1 IS. 


UBV I 


J 06 


9 no of a 


UBVRI 


J 06 


2001V 


UBVR 


V03 T,06a 


2001ba 


BVI 


K04 


2001bt 


BVRI 


K04 


2001cn 


BVRI 


K04 


2001cz 


BVRI 


K04 


2001el 


BVRI 


K03 


2002bf 


BVRI 


L05 


2002bo 


BVRI 


K04 


2002or 


UBVR 


P04 


2003cg 


UBVRI 


ER06 


2004S 


UBVRI 


K07 


2004eo 


UBVRI 


P07b 


2004fu 


BV 


T06 


2005am 


BVR 


L06a 


2005cf 


UBVRI 


P07a 


a 


Photometry 


' References: 



H91 IHamuv et all (1991), 
P87: IPhillips etal] 11987), 

W94: IWells et al l jHm ), 
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of the model directly when applying the model to SN 
observations. Example fits to high and low redshift SNLS 
SNe are shown in figures [5] and [5] (the same SN were used 
as examples in GO 7). 

The parameters of the color relations used to determine 
parameter C (equation[2|) are given in tables [3] and [U and 
the U — B and J/02 — B relations shown in figures ITD1 and 
[TT1 The U - B relation is derived from 72 SNLS and 33 
low-z SNe, and the U 02 - B relation from 63 SNLS SNe. 
Here, because there are few SNe with both V and J/02 
measurements, we include data with no rest-frame V ob- 
servations by using their U — B colors to predict B — V, 
incorporating the additional uncertainty due to o~ c \ T and 
the appropriate covariances. The requirements to be in- 
cluded in the color relation fit are less stringent than 
those required to be used when measuring 0o since we 
do not require that the fits to each filter be strongly con- 
strained in isolation. The primary limitation is instead 
that most SNe do not have simultaneous measurements 
oiU-B and B — V. Carrying out this test with SALT2 
gives similar values for the cr c i r parameters. It is interest- 
ing to consider the residuals from the U — B versus B — V 
relation. The slope of the dominant relation is clearly in- 
consistent with that expected from extinction. However, 
there are a number of outliers that are displaced in a di- 
rection that could be explained by extinction. This leads 
us to hypothesize that the color relations for the majority 
of SNe in our sample are dominated by intrinsic effects, 
with a minority instead dominated by dust. 

Following G07, we can get some feel for how useful 
extending the model blue-ward of U is be comparing 
measurements of C if the model is cut at different wave- 
lengths. This is shown in table [5] for 04D3cy (z — 0.64) 
and 04D3gx (z — 0.91). The errors include statistical 
uncertainties, template errors, and X-correction errors 
as described in H07. The benefits for the higher redshift 
SN are clear. 

The stretch is only defined relative to some fiducial 
template, and does not have any independent meaning 
when comparing different fitters, or even different ver- 
sions of the same fitter. Instead, a variable like Amis 
can be used which has a more physical definition. 8 How- 
ever, in practice, because light-curve fits are relatively 
insensitive to the exact shape of the light-curve template 
near peak, even Amis is somewhat ambiguous, and will 
change when the light-curve model is updated. There- 
fore, we simply provide the relation between Amis and 
sb for our template: 

Amis = 1.00-1.63 (s B - l)+2.03 (s - 1) 2 -1.82 (s - l) 3 . 

This relation should not be used for precision work, such 
as attempting to combine data fit with different packages. 

5.1. Comparing the color parameters from SNLS and 
low-z data 

For the U — B color relation, we can split the sam- 
ple into low-z and SNLS data as a test of SN evolu- 
tion. This results in the values given in table [3J There 
is some tension between the values. However, this only 
includes statistical errors, and systematics are expected 

8 The quantity Amis is the amount in magnitudes which the B 
band light curve has dimmed 15 rest frame days after peak. 



TABLE 3 
SiFTO Color Relation 
Parameters 



Parameter 


Value Error 




SNLS 


± low-z 



a 0.352 0.049 

b 0.218 0.064 

c -0.051 0.007 

cr cll 0.055 0.005 



SNLS only 

a 0.435 0.067 

b 0.302 0.076 

c -0.067 0.008 

o- clr 0.044 0.006 



low-z only 

a 0.259 0.090 

b 0.145 0.132 

c -0.037 0.016 

a ch 0.071 0.011 



Note. - The coefficients 



for the U — B vs. B — V color 
relation of equation [2] Only 
the stati stical errors are given. 
See i]5.1l for a comparison of the 
SNLS and low-z values. 



TABLE 4 
SiFTO Color Relation 
Parameters: U02 — B 



Parameter Value Error 



a 02 0.121 0.050 

602 0.172 0.138 

c 02 -0.110 0.027 

°02cir 0.050 0.008 



Note. - The coefficients 



for the U02 — B vs. B — V color 
relation. 



TABLE 5 

Effect of Including Near-UV Information 





bands 


C 


WU02B a 


WUB b 


WBV C 








SNLS-05D31b 2 = 


= 0.65 






3980 A 
3250 A 
2700 A 


IMZM 

tm^mzm 

gMTMi-MZM 


0.066 ± 0.097 
0.008 ± 0.040 
-0.002 ± 0.038 





0.43 




0.69 
0.40 


1 

0.31 
0.18 






SNLS-04D3gx 2 : 


= 0.91 







3980 A i M z M -0.106 ±0.080 1 

3250 A r M iMZM -0.100 ±0.061 10 
2700 A g M r M iMZM -0.097 ± 0.040 0.60 0.40 0.00 



Note. — The er 


rors in C for two SNLS SNe when the model is cut at different W 




lengths. The cuts • 


correspond to including all of Uq 2 U BV , removing U Q2 , and remc 


ving 


U. Note that the rc 


ist-frame V is not measured for SNLS-04D3gx. The improvement ii 


1 the 


measured color for 


04D3gx when r is added is due to the improved constraints or 


1 the 


other light-curve p, 


iramcters such as the date of maximum, even though it is not forn 


rally 


used in the color. 






a The weight given 


to UQ2 ~ B when forming C. b The weight given to U — B when for: 


ning 


C. C The weight giv 


en to B — V when forming C. 
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to be of similar importance. The details of our system- 
atics analysis will be presented with our cosmological re- 
sults elsewhere, but it is useful to summarize some of 
the more important contributions. The most important 
uncertainties arise from the uncertainty in the magni- 
tudes of our flux standard (Vega) on the Landolt system 
(which affect our if-corrections) and the uncertainty in 
the calibration of the low-z data. The former are dom- 
inate d by the raw observation al uncertainties of Vega 
from I Johnson fc Morgan! (|1953| ). To determine how well 
calibrated the low-z data is to the Landolt system, we 
compare photometry of the same SNe obtained by dif - 
ferent observers (e.g., iRiess et alJll999HJha et al.ll2006fl , 
and conclude that there is a global zeropoint uncertainty 
of about 0.015 mag in BV for each survey. For U the situ- 
ation is more complicated. As noted previously, observer 
frame U data is more difficult to calibrate than BV, so 
we double this value to 0.03 mag in this filter. However, 
as it turns out the B zeropoint uncertainty dominates 
over U, at least partially because it affects both U — B 
and B — V . 

In order to test the effects of these systematics, we ad- 
just the relevant parameter (e.g. the U zeropoint), then 
re-derive the SiFTO model and color relations. We then 
consider the change in osnls — aiow-z to be the resulting 
systematic in a, etc. For the uncertainties in the magni- 
tude of Vega, we find a a = 0.014, <Jb = 0.050, and a c = 
0.007. We do not include cr c i r because its value is likely 
to simply test the homogeneity of the data sample rather 
than be a good test for evolutionary effects. For the zero- 
points of the low-z data, we find a a = 0.022, erf, = 0.022, 
and cr c = 0.013. Other, smaller, contributions include 
the uncertainty in the MegaCam and Landolt bandpasses 
and the MegaCam zeropoints. The latter is particularly 
important for b, giving ct& = 0.034 for sub-1% zeropoint 
uncertainties. 

Putting this information together, we conclude that 
5a = 0.179 ± 0.112(stat) ± 0.030(sys), 5b = 0.157 ± 
0.152(stat) ± 0.068(sys), and 5c = 0.030 ± 0.018(stat) ± 
0.019(sys). The resulting x 2 is 4.57 for 3 degrees of free- 
dom, which is expected to occur by random chance in 
the case of no evolution 21% of the time. 

6. HOW CLOSE TO PEAK DO WE NEED DATA? 

One question we can ask of our model is whether or 
not it can predict the peak luminosity of the SN when 
we have no observations at that epoch. That is, how well 
can our model extrapolate from late time observations to 
the peak? In principle, the template error snake should 
allow us to handle situations where there is little or no 
data near peak, since the result will simply be a measure- 
ment with very large errors. However, the error snake is 
not known perfectly, so we would like to test this state- 
ment. Most SN cosmology analyses have required that 
SNe have data near peak, but the exact value used for 
this cut varies from paper to paper and is usually not 
explained. The previous generation of SN fitting pack- 
ages, which often did not incorporate template errors, 
frequently produce implausibly precise results when used 
in this situation. 

Our test is based on real observations; a test on sim- 
ulated data is currently likely only to be useful as an 
internal consistency check. We collect a sample of very 
well sampled, high-S/N light-curves with data well be- 



fore peak, and then, for each SN, compare the results 
of the fits to all of the available data to those where we 
steadily remove more and more photometry from around 
the peak. We parameterize this by the daygap, defined as 
the epoch of the first included point relative to the date 
of B maximum, and explore the range from —5 to +15 
days. Our sample consists of the 11 low-z SNe: 1992 A, 
1995D, 1999aa, 1999dq, 1999ee, 2000E, 2002bo, 2002er, 
2003du, 2004eo, and 2005cf. Even these SNe generally do 
not have data before day —10, which is why we have not 
tried to answer the complementary question of how well 
our model can predict the peak magnitudes of SNe la if 
we only have pre-maximum observations. 

There are two things we want to test: that the results 
are unbiased, and that the errors reflect the missing data 
appropriately. For the latter, we compare the x 2 °f the 
difference between the light-curve parameters for the full 
light-curve, and for those at various values of the daygap. 
Since the different fits to the same SN should be posi- 
tively correlated, we will probably slightly over-estimate 
our errors. This test is sensitive to small calibration 
difficulties that arise when different telescope/detector 
systems observe different parts of the light-curve, which 
is frequently the case with low-z data. This will artifi- 
cially inflate the differences. Therefore, the precise value 
should not be taken too seriously; what we are looking 
for are sharp changes in the x 2 - 

The most important parameter to test de- 
fined in equation[TJ which is used to construct the Hubble 
diagram. Rough numerical simulation suggests that bias 
values larger than about 0.02 are significant. Most of 
the variables show relatively little structure in their x 2 , 
so we only present the two given in figure [T2l First, the 
X 2 of the peak magnitude in B rises sharply around day 
5-10 after peak. Second, while this effect is also visible 
for m corr , the effects are slightly muted. Our tests show 
little evidence for bias in the corrected peak magnitude. 
Based on the x 2 -hke test, a good policy is to require data 
within ~ 7 days of peak when using SiFTO. 

Note that the accuracy of the derived parameters does 
degrade considerably as the daygap is increased in all 
parameters (see the bottom of figure [T2]) . so SNe with 
no data near peak will contribute relatively little infor- 
mation to the Hubble diagram. We carry out the same 
test with other fitters in H7.21 and find generally similar 
results. 

7. COMPARISONS WITH OTHER FITTERS 

In this section, we compare the results of SiFTO with 
those of other fitters, particularly the published versions 
of SALT and SALT2. Our main focus will be on compar- 
isons of SiFTO, SALT2, etc. when they are used as dis- 
tance estimators rather than as light-curve fitters. There 
are a number of non-SN related settings that affect this 
comparison (i.e., the assumed filter response of the Lan- 
dolt filters 9 , etc.), so this comparison is not entirely fair. 
In § E] we present the results when these differences are 
removed which provides a much more accurate test of 
how different the results of the fitters really are. 

We discuss MLCS2k2 in a more limited fashion in this 
section, keeping in mind the following limitations: First, 

9 SALT a nd SA LT2 follow A06 in using blue-shifted versions of 
the [5 esse 1 n irr99oh responses, while SiFTO does not. 
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MLCS2k2 makes use of a number of priors which should 
be accounted for in any detailed comparison, and which 
in principle should be modified for the particular data set 
under consideration. Second, it fits in magnitude space, 
which will introduce a bias when comparing low signal- 
to-noise data 10 . Finally, the current version of MLCS2k2 
does not properly handle the case where more than one 
observed filter maps to the same rest frame Landolt filter. 
This limits the sample size that can be compared using 
SNLS data. 

We also consider comparisons of the Hubble diagram. 
All of the fitters require that some additional uncertainty 
be added to the luminosity distances when computing 
the best-fit cosmology. This is often referred to as "in- 
trinsic scatter" or "intrinsic dispersion," and presumably 
reflects our imperfect understanding of SN la. We de- 
note this quantity by cr; n t, and find a mt = 0.11 mag for 
SiFTO. Comparing the values of <7; n t for the different fit- 
ters is not currently a useful mechanism for comparison. 
The issue lies with the limitations of the template error 
snakes, which are difficult to determine and are very sen- 
sitive to how accurate the input flux errors are measured. 
In addition, currently both SiFTO and SALT2 do not 
include off-diagonal terms in the error snake (i.e., differ- 
ent epochs are considered uncorrelated), an assumption 
which is clearly not true. This is particularly an issue for 
SALT2, which imposes the color relation during the fits. 
Both effectively absorb this unknown into <7i n t , and there 
is no reason to expect this to affect both equally. There- 
fore we limit ourselves to the more robust comparison of 
the RMS of the resulting Hubble diagram. However, the 
RMS also includes the effect of measurement errors, so 
is only useful as a comparison statistic if precisely the 
same set of input data are used for all fitters. A com- 
plication is that there is considerable overlap between 
the data used to measure the RMS and that used to de- 
rive the SiFTO model; the same issue applies to SALT2, 
MLCS2k2, and (to a lesser extent) SALT. However, we 
note that SiFTO /SALT / S ALT2 make absolutely no use 
of relative distance information as part of the training 
process. 

7.1. Comparison of SiFTO and SALT2 light-curve fits 

We begin by comparing SALT2 and SiFTO as light- 
curve fitters - that is, in terms of how well they fit the 
available photometry rather than in terms of the pre- 
cision of their relative distances. Here we do have the 
problem that the training process is specifically designed 
to reduce the scatter of the input photometry with re- 
spect to the model, so it is beneficial to train and test 
the models with disjoint data sets. We have carried this 
process out for SiFTO by randomly splitting the training 
sample in half (giving about 90 SNe for both testing and 
training), training a set of templates on the first half, 
and then applying that model to the second half. We 
then look at the residual between the model and input 
flux normalized so that the predicted peak flux in each 
filter is unity. Because this includes observational errors, 
we split this up into redshift bins (Az = 0.1). In order to 
mitigate photometry failures, we sigma-clip the residuals 

10 Subsequent to this writing, a version of MLCS2k2 which works 
in flux space was kindly provided by G. Narayan (2007, private 
communication). However, it does not address the other issues. 



TABLE 6 
Comparison of Low-z Hubble 
Diagram Fits 



Method 


RMS (mag) 


o-int (mag) 


SALT 


0.156 


0.134 


SALT2 


0.158 


0.129 


SiFTO 


0.155 


0.119 



NOTE. — The RMS and cr int values for fits to 
the same low-z data, and assuming ilrn, = 0.3, 
f2^Y — 0.7. ^ijxt * s n °t useful for evaluating the 
different fitters. 



at 4c for both fitters. The results are shown for SALT2 
and SiFTO in figure [13] using the same photometry. This 
comparison is unfair in two ways: first, the splitting into 
disjoint training and testing samples has not been car- 
ried out for SALT2. Second, the SiFTO model has more 
free parameters since a color relation is not imposed dur- 
ing the fitting process. Nonetheless, the two fitters give 
comparable results. 

7.2. Comparisons using low-z data 

The different light-curve fitters generally agree quite 
well on the derived light-curve parameters, particularly 
when the MCLS2k2 priors are removed. Examples are 
shown in figures [TJ] and [T5l We generally do not expect 
perfect agreement between the different light-curve pa- 
rameters, but only that there be a clear relation between 
the values. This is particularly true when discussing mea- 
sures of the light-curve shape, which are defined very dif- 
ferently, and SN color, which can have an arbitrary offset 
related to the intrinsic color of the fiducial SN la. The 
data set for this comparison is 78 low-z SNe from the 
literatu re, using the same general requirements as were 
used in lConlev et al.l (|2007f ). although the sample here is 
larger because we can include SNe which are not in the 
smooth Hubble flow. 

The derived parameters are clearly quite correlated be- 
tween the different fitters, although there are a number 
of outliers in each relation. The ms comparison is con- 
sistent at about the 1% level globally, and since for most 
purposes only the relative distances matter, this is en- 
couraging. 

Next we turn to the RMS around the low-redshift Hub- 
ble diagram. We take the above sample, remove all SNe 
with z < 0.015, and then fit the Hubble diagram as- 
suming fl m = 0.3 and Qa — 0.7 - the exact values are 
unimportant because of the limited redshift range of our 
test sample. The RMS values and <7i n t are given in ta- 
ble [51 MLCS2k2 is not included in this test because of 
complications related to the interpretation of SN colors, 
which in turn is related to the question of whether or not 
the low-z Hubble flow is really smooth or has a Hubble 
bubble, somet hing that MLCS2k 2 favors and the other 
fitters do not (jConlev et al.ll2007h . Since we do not in- 
clude a bubble in our test, this will tend to artificially 
inflate the RMS for MCLS2k2. However, we do perform 
the test using a larger data set below, where this issue is 
mitigated by a larger, higher-redshift sample. 

We find that the three fitters perform equally well on 
this data set. This is simultaneously encouraging and dis- 
couraging; it gives us confidence in the different fitters, 
but at the same time suggests that it is not easy to signif- 
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icantly improve the quality of the distance estimates for 
SN la beyond the 7 — 9% level without some improve- 
ment in the training data set. Interestingly, when we 
perform the same comparison with high-z SNe, SALT2 
and SiFTO do considerably better than SALT. 

Finally, we can test the other fitters to see how well 
they fare if there are no data near peak, using the same 
method and set of 11 SNe as in § [6l It is particularly 
interesting to consider SALT, which alone among those 
discussed here does not have a template error snake. The 
bias values are shown in figure [16] and the x 2 -like test 
results in figure H7l MLCS2k2 performs the best in the 
X 2 -like test, which is also true for the individual fit quan- 
tities. MLCS2k2 and SALT display some bias even for 
low values of the daygap. One possible explanation is 
that both fitters would benefit from some improvement 
in their pre-maximum light-curve model. 

We note that SALT does quite poorly in the x 2 -like 
test, as expected. Furthermore, we obtain similar results 
when we do not use the template error snake in SiFTO. 
The real lesson is then that including template errors 
plays a critical role in obtaining reliable fits when there 
are no data near peak. 

7.3. Comparisons using high-z data 

It is also interesting to compare the results of fits to 
high-redshift SNLS data. We first compare SiFTO to 
the published versions of SALT and SALT2. Because 
we can apply each fitter to the same photometry, again 
the best comparison statistic is the RMS of the result- 
ing Hubble diagram around the best fit. The alterna- 
tives all have problems. In particular, comparing the re- 
sulting errors on the cosmological parameters is a badly 
flawed metric unless the fits are virtually identical, given 
the non-linearity of the parameter space. An additional 
complication occurs if the constraint contours impinge 
on a physically prohibited region, such as Q m < 0. Tests 
with synthetic data sets are perhaps useful as an inter- 
nal consistency check, but we do not currently under- 
stand SNe la well enough for them to have any absolute 
meaning. 

For this set we start with the set used in A06 and mod- 
ify it slightly. We removed two SNe without solid spec- 
troscopic type identification (04D2iu and 04D3is), one 
with photometry problems (03D4cn), and one which is 
probably SN 1991bg-like, and hence not well represented 
by any of the models (SN 1992bf). We also removed a 
number of SNe whose light-curves are poorly sampled 
(SNe 1997Y, 1998eg, 03D3bh, 03D4gl, 04Dlag). Finally, 
we added two low-z objects, SNe 2001V and 2002bf. This 
leaves us with a sample of 43 nearby SNe and 65 high- 
redshift ones. 

We carry out three cosmological fits for each: a flat- 
Universe Q m fit, a Q m + Q\ fit, and a flat-Universe 
w fit. For the latter, we incorporate the SDSS 
baryon acoustic oscill ations (BAO) measurements of 
lEisenstein et al.l ([2005). since these are marvelously or- 
thogonal to the SN constraints. The results are given in 
table [71 A bootstrap analysis indicates that the error in 
the RMS for this sample is ~ 0.005 mag. From the per- 
spective of the RMS test, SiFTO and SALT2 are indis- 
tinguishable, and both considerably outperform SALT. 
We also demonstrate the effects of not incorporating the 
U02 filter in SiFTO, which worsens the RMS values at 



TABLE 7 

Comparison of Fits to A06 Photometry 



Fit type RMS Parameters 



SiFTO 



Flat £lrn 


0.161 


°2m 


= 0.263 ±0.034 


$i m -\- Q A 


0.160 




+ Q A = 1.358 ±0.353 








-n A = -0.550 ±0.108 


Flat f2 m + w 


0.161 




= 0.274 ±0.021 






W = 


-1.047 ±0.088 



SiFTO w/out U 2 



Flat U m 




0.171 


— 


0.261 ±0.035 






0.170 


± 


Q A = 1.438 ±0.366 








- ' rri 


Q A = -0.568 ±0.110 


Flat O m + 


w 


0.172 


\lm — 


0.273 ±0.022 








UU = - 


-1.054 ±0.088 


SALT 


Flat dm 




0.187 


ilrn — 


0.260 ±0.042 






0.187 


± 


Q A = 1.023 ±0.388 










C A = -0.477 ±0.120 


Flat Q m -f 


w 


0.188 




0.272 ±0.022 










-1.030 ±0.098 


SALT2 


Flat n m 




0.160 


' ' — 


0.222 ±0.034 






0.160 


± 


n A = 1.306 ±0.350 










Q A = -0.613 ±0.100 


Flat H m + 


w 


0.160 


\lm — 


0.263 ±0.021 








W = - 


-1.126 ±0.090 



NOTE. — RMS values (in magnitudes) and resulting cosmo- 
logical parameters for fits to the A06 photometry using SiFTO, 
SALT, and SALT2. Only statistical errors arc included. The 
results of applying the SiFTO model limited to A > 3300A (ex- 
cluding c^Q2) arc also shown. As explained in the text, only the 
RMS should be used as a performance comparison statistic. The 
RMS has an uncertainty of ~ 0.005 mag. 



approximately the 2a level. We note that the benefits 
of extending the model bluer than U are greater for the 
third year SNLS sample than are shown here. This is 
probably a result of the fact that the zm exposure times 
have been increased, so our z > 0.8 SNe now have better 
measurements of tub and hence have increased weight in 
our fits. The RMS including high-z data is slightly worse 
than if only low-z SNe are used, which simply reflects the 
larger observational errors of the distant sample. 

Some of the resulting constraints are shown in fig- 
ure [THJ We see that some differences remain between 
the fitters. We shall turn to the question of whether this 
reflects significant differences between the models in fJSJ 
A weakness of the RMS as a test statistic is that it is sen- 
sitive to outliers. Including 3er iterative outlier rejection 
does not appreciably change the results. 

We can perform a more limited version of this test in- 
cluding MLCS2k2, although given the caveats given ear- 
lier this is somewhat less certain and has a smaller sam- 
ple size. We obtain t he results of MLCS2k2 fi ts to the 
A06 photometry from lWood-Vasev et all (|2007f) . This is 
a subset of the full data set, with about 93 SNe in total. 
We therefore limit the SALT, SALT2, and SiFTO fits to 
the same sample. We do not include the ESSENCE SNe 
in this test. Carrying out a similar analysis as the above, 
the RMS for SALT is 0.179 mag, for SALT2 0.159 mag, 
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for SiFTO 0.160 mag, and for MLCS2k2 0.205 mag. As 
noted for the low-z comparison, the MLCS2k2 value is 
somewhat inflated by the issue of the Hubble bubble. If 
we approximately remove this effect by hand, the RMS 
for MLCS2k2 is about 0.19 mag, so we conclude that its 
RMS is similar to SALT. Adjusting the priors to better 
match the survey properties of SNLS might improve the 
performance of MLCS2k2. It will be interesting to per- 
form this comparison again once the issues surrounding 
SN colors have been resolved. 

8. A BETTER COMPARISON BETWEEN SIFTO 
AND SALT2 

The differences between the SALT2 and SiFTO results 
on the A06 photometry reflect not only differences in the 
underlying model but also a variety of incidental settings, 
many of which are not related to SNe at all. While these 
differences are interesting, and play a role in the system- 
atics error budget of our final result, here they simply 
obscure what we really want to know: what effect do the 
differences in approach have on the cosmological param- 
eters? 

We can try to answer this question better by using 
the same non-SN related settings for both fitters (filter 
functions, the magnitudes of Vega on the Landolt sys- 
tem, etc.). The published version of SALT2 was trained 
on the set of SNLS photometry available at the time us- 
ing the MegaCam calibration of A06. There are now 
more SNLS SNe available, and the calibration has been 
improved. We can further enhance the comparison by 
training both fitters on exactly the same low-z and SNLS 
photometry, using the same filter functions assumed by 
SALT2 in SiFTO, and adding the additional low-z SNe 
provided by T. Matheson to SALT2. 

Carrying out this prescription, we obtain the results in 
table [U and figure [T9l The results substantially closer to 
each other. This is very encouraging, and suggests that 
in terms of the cosmological results the two approaches 
agree very well. 

9. CONCLUSIONS 

We have described SiFTO, a light-curve analysis pack- 
age developed for use with SNLS data. It is most similar 
in spirit to SALT/SALT2: all three work by adjusting 
a model of the SED rather than with templates in some 
rest-frame filters, and, like SALT2, SiFTO includes high- 
z SNLS data in the training process. This allows us to 
extend our model blue-ward of the rest-frame U, which 
is extremely valuable when fitting high-z SNe. Distances 
and residuals from the Hubble relation are not used in 
our training process, so this is not dependent on the cos- 
mological parameters. The model can be applied to ob- 
servations from 2700 to 7100 A in the rest frame. 

The greatest difference between SALT/SALT2 and 
SiFTO is the manner in which SN colors are handled. 
SiFTO does not impose a color model during the actual 
light-curve fit, but instead adjusts the SED to match the 
observer frame colors. The SED is then used to mea- 
sure various rest-frame colors, and these are combined 
to form a single color parameter C. The relation used 
to combine measurements is derived using SNe at inter- 
mediate and low redshifts where each of the colors are 
well measured. The advantages of extending the model 
further towards the near-UV, made possible by including 



TABLE 8 

Comparison of Fits to A06 Photometry with 
Retrained Models 



Fit type RMS Parameters 



SiFTO 



Flat H m 


0.159 


— 


0.262 ±0.033 




0.161 


Sz m ■+■ 


C A = 1.275 ±0.357 








Q A = -0.536 ±0.106 


Flat £~2 m + w 


0.160 


O m = 


0.272 ±0.021 






w = - 


-1.042 ±0.086 


SALT2 


Flat £7 m 


0.159 


— 


0.254 ±0.034 




0.160 




Q A = 1.393 ±0.353 








Q A = -0.578 ±0.112 


Flat £~2 m + w 


0.160 


Vim — 


0.271 ±0.021 






W = - 


-1.064 ±0.088 


Note. — R.MS 


values (ir 




ides) and resulting cosmo- 


logical parameters 


for fits to 


the A06 


photomctry using versions 


of SiFTO and SALT2 train* 


id on the 


same data with identical 


non-SN settings. The SiFTO 


values ai 


•e slightly different than in 


tablc^jbccause her 


e SiFTO il 


5 using th 


e MegaCam filter responses 


used by SALT2. 









hi-z data, are significant. This is an argument in favor 
of the light-curve fitter class of packages (such as SALT, 
SALT2, and SiFTO) versus distance estimators, which 
cannot incorporate high-redshift data as easily. 

We find that SiFTO produces the most reliable results 
when photometry is available within 7 rest frame days of 
the B peak, and data prior to peak is not necessary. The 
same result applies to the other fitters as well. Including 
a template error in the model is particularly important 
when there are no data near peak. 

We have carried out a number of tests comparing 
SiFTO to other packages. We find that SALT, SALT2, 
SiFTO, and MLCS2k2 agree fairly well in terms of the 
derived light-curve parameters when the full SN sam- 
ple is considered. In terms of their performance on 
the Hubble diagram, for low-z data SALT, SALT2, and 
SiFTO are comparable. However, when the test is 
extended to higher redshifts, SiFTO and SALT2 out- 
perform SALT, and SiFTO performs as well as SALT2, 
while MLCS2k2 performs about as well as SALT. These 
comparisons show that there are real, although small, 
differences in the derived cosmological parameters from 
each when they applied to the same photometry. If we 
train SiFTO and SALT2 using the same non-SN related 
settings and photometry, we find that the differences are 
much smaller. This indicates that, for currently available 
data, the fundamental differences between the SALT2 
and SiFTO models are minor despite the significant dif- 
ferences in approach and algorithmic design. 
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06D2ag, z=0.310, sifto, spline 




3000 4000 5000 , 6000 

Wavelength (A) 



Fig. 1. — Example of adjusting the SED to have the desired observer frame colors using the relatively low-z SNLS SN 06D2ag (z = 0.310). 
Here Zm ls n °t used because it is redder than the range over which the model is trained. The dotted line is the fiducial template tj>o, the 
solid red line is the adjusted SED, the smooth line is the warping spline, the open circles are the original integrated flux through each of 
the observed filters (also shown), and the filled circles are the final, target fluxes. The rest-frame magnitudes are evaluated by performing 
synthetic photometry on the SED. The wavelengths are given in the rest frame. 
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Fig. 2. — Difference between sb and Sj as a function of sg for all the SNe in the training sample. SNLS SNe la are shown as black 
circles, low-z SNe la as blue circles. The best fitting linear model (including errors in both sb and Sj) for both data sets is over-plotted. 
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Fig. 4. — Stretch in a particular observed filter (sj) as a function of sg for various effective observed frame wavelengths Xj . The MegaCam 
filter %m is a good match to B at z = 0.75, U at z = 1, and is between B and V at z = 0.5. 
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Fig. 5. — Filters used when constructing the SED template ipo (solid lines) compared with UBVRI (dashed lines). The B and / filters 
(dot-dashed lines) are shared by both. The normalization of the filter responses is arbitrary. 
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Fig. 6. — Example error snakes for some of the template filters. ARB1 is bluer than U and L0G7 is the bluer of the two filters that 
overlap with R. Both ARB1 and LOG7 are synthetic filters. 
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Fig. 7. — Light-curve template fits in our filters. The red squares are SNLS data, and the blue circles are low-z data. The fluxes have 
been normalized so that the B template has a peak value of unity, and the colors to be representative of a typical SN la. A full resolution 
version is available at http://qold.astro.utoronto.ca/conley/sifto_fullres.pdf. 
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Fig. 8. — Sample SiFTO fit to 04D3gx, a high-redshift (z = 0.91) SNLS SN la. The dashed lines are the error snake for the template, 
and do not include the uncertainties in the light-curve parameters. This SN is from the first-year sample, and therefore does not represent 
the quality of light-curves currently obtained by SNLS. In particular, the signal-to-noise ratio in zm has improved considerably in the 
third-year sample due to longer exposure times. 
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Fig. 9. — Sample SiFTO fit to 03D4ag, a low-redshift (z = 0.285) SNLS SNe la. The dashed lines are the error snake for the template, 
and do not include the uncertainties in the light-curve parameters. 
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Fig. 10. — Derived U — B vs. B — V relation using low-z (blue circles) and SNLS data (red squares). The best fit relation to the full 
sample is shown as a solid line. Also shown is the relation one would expect from Milky- Way like dust with Ry = [1.6, 3.1, 4.6]. The points 
rejected as outliers are open. 
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Fig. 12. — Reduced \ 2 f° r th e corrected magnitude (equation [TJ as a function of the gap between the first data point and the date of peak 
B luminosity (top panel). The middle panels show the same relation for the raw B magnitude and the bias in the corrected magnitude, 
and the bottom panel shows the error in the B magnitude. 
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Fig. 13. — Comparison of the RMS around the light-curve fits as a function of redshift for SALT2 (red squares) and SiFTO (blue 
circles). The residuals are normalized so that the model peak flux in each filter is unity. Note that the SALT2 training sample includes the 
photometry used to perform this test, and that SiFTO has more free parameters in the fits. 




Fig. 14. — Comparisons of the light curve shape parameters for SiFTO vs. SALT2 (left) and MLCS2k2 (right). Note that the meanings 
of these parameters are quite different. Linear fits are given by black dashed lines, and 3 rd order polynomial fits as red dot-dashed lines. 
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Fig. 15. — Comparisons of estimated peak colors for SiFTO vs. SALT2 (left) and MLCS2k2 (right). The definitions of the peak color 
generally differ by a constant offset, so here we are testing if they are linearly related with a slope of one. The dashed lines are the best fit 
linear relations, which have slopes of 1.04 ± 0.02 and 1.2 ± 0.03, respectively. 
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Fig. 16. — Bias in the corrected peak magnitude relative to the case when all available data are used vs. the daygap for SALT, SALT2, 
and SiFTO, as well as for the equivalent MLCS2k2 variable (/i). The error in this measurement is ~ 0.015 mag. 
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Fig. 17. — Reduced x 2 f° r the corrected peak magnitude m CO rr (i.e., the value used to construct the Hubble diagram) vs. the daygap 
for SALT, SALT2, and SiFTO. Also plotted is the reduced \ 2 f°r the MLCS2k2 distance modulus fi, which is roughly equivalent. 




Fig. 18. — Comparison of SiFTO, SALT and SALT2 fits to a subset of the A06 photometry. Only the 68.3% confidence limits are shown 
for clarity. In addition, the constraints the result when (7o2 is not included are also shown. 
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Fig. 19. — Comparison of SALT2 and SiFTO with the same non-SN ancillary settings and trai ned on the same photometry. This is a 
better representation of the fundamental differences between the models than that shown in figure 1181 



